Realization of quantum gates with multiple control qubits or multiple target qubits in 

a cavity 
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In this paper, we propose a scheme to realize three-qubit controlled phase gate and multiqubit 
controlIed-NOT gate of one qubit simultaneously controlling n target qubit with four level quantum 
system in a cavity. Adjustment of level spacing during the gates implementation is not required. 
Implementation time for multiqubit controlled NOT gate is independent of number of qubit. Three- 
qubit phase gate is generalized to n-qubit phase gate with multiple control qubit and the number 
of steps (complexity) reduces linearly as compare to conventional gate decomposition method. Our 
scheme can be applied to various types of physical systems such as superconducting qubits coupled to 
resonator and trapped atoms in a cavity. Experimental possibility of our approach is also presented. 
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I. INTRODUCTION 



Quantum computing has the potential ability to carry 
out certain computations much faster than classical com- 
puting. For example factorization of a large number 
via Shor's algorithm [l| and the search of an item in 
an unsorted database containing N elements Two- 
qubit gates and one-qubit gates are the building blocks 
for quantum computing networks [^. Many physical 
systems have been proposed as candidates for imple- 
mentation of quantum information processing like atom 
in cavity quantum electrodynamics (QED) and nuclear 
magnetic resonance (NMR). Among them cavity QED 
analogs with superconducting qubit systems are getting 
favorable attention [3| . A two-qubit gate was experimen- 
tally realized using superconducting qubit systems cou- 
pled through capacitors 0, 0] , mutual inductance , or 
cavities M- 



Multiqubit quantum gates play a significance role 
in quantum information processing. Experimentally, a 
three-qubit controlled NOT gate has been demonstrated 
with trapped ions @. Multiqubits gates constructed 
through conventional gate decomposition method (lo| . 
usually makes the procedure complicated for the case of 
large number of qubits (for example see Sec. III. A). This 
is because the number of single qubit gate and two qubit 
gates required for the implementation heavily depends 
on number n of the qubit. The purpose of this work is to 
realized three-qubit controlled phase gate and multiqubit 
controlled NOT-gatc of one qubit simultaneously control- 
ling n qubits (which we will denote as NTCNOT-gate) 
in a cavity QED. Lets first introduce these gate below 



A. Two kind of multiqubit quantum gates 

In three-qubit quantum controlled phase gates when 
two control qubits \qi) and \q2) are in state |1), phase 
shift e^^ induces to the state |1) of the target qubit jga). 
When control qubits are in state |0) nothing happens to 
target qubit . This transformation can be written as [llj | 
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Here, Sq-^,i, Sq^.i, and Sq^^i are the standard Kronecker 
delta functions and \qi) , \q2) and jga) stand for basis 
states |0) or |1) for qubit 1, 2 and 3. Circuit for three- 
qubit controlled phase gate is shown in Fig. Ufa). Thus 
three-qubit quantum phase gate introduces a phase 77 
only when the input state of all three qubits in |1) . In 
this proposal, we discuss the implementation of three- 
qubit quantum phase gate with rj = n. It may be men- 
tioned that three-qubit controUed-NOT gate (known as 
Toffoli gate) can also be achieved using present proposal. 
Toffoli gate is equivalent to three-qubit controlled phase 
gate plus two Hadamard gates on target qubit as shown 
in Fig. mh). 

Here, we consider NTCNOT-gate consist of control 
qubit 1 and n target qubits labeled by 2,3, n as shown 
in Fig. Ilja). We define control qubit in |0), |1) basis and 
each target qubit in |+), |— ) basis. Thus, input state can 
be written as 



10) IK 

k=2 



c)+ii)n( 



>fe)- (2) 



k=2 



When NTCNOT-gate is applied to above state, we ob- 
tained final state given by 
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FIG. 1: (a) Three-qubit controlled phase gate. Z represents 
Pauli rotation CTz. If control qubits 1 and 2 (shown by filled 
circle) are in state |1) then only state |1) at Z is flipped by 
phase of vr. When control qubits are in state |0) nothing 
happens to target qubit. (b) Relationship between three- 
qubit CNOT-gate (known as Toffoli gate) and three-qubit 
controlled phase gate. The circuits on left side and right side 
of (b) are equivalent to eachother. The symbol ® on left side 
of (b) represents NOT gate on target qubit. If control qubits 
1 and 2 are in state |1) then state at ® is flipped as |1) — ^ |0) 
and |0) — >■ However, when control qubits 1 and 2 are in 
state |0) then state at © remain unchanged. For right side of 
(b) portion inclosed in dashed box represents three-qubit con- 
trolled phase gate. The element H is called Hadamard gate 
and leads the transformation |0) --^ ]+) = (1/V2){\0) + 
and|l)^|-) = (l/y2)(|0)-|l)). 



It is clear from Eq. [2] and |3] that when control qubit is 
in state |1) then state at each target qubit is flipped as 
|-|-) — )■ |— ) and |— ) — >■ |-f) . If control qubit is in state |0) 
nothing happened to each target qubit. 

It may be mentioned that NTCNOT-gate can be de- 
fined in |-|-), |— ) basis. However, two Hadamard gates 
on control qubit before and after the phase gate of one 
qubit simultaneously controlling n target qubits would 
be required as shown in Fig. [2](b). The NTCNOT-gate 
can also be defined in |0), |1) basis. However, in this case 
Hadamard gate on each target qubit before and after the 
NTCP-gate (i.e., total 2(n — 1) Hadamard gate) would 
be required as shown in Fig. [2](c). In contrast, defining 
control qubit in |0), |1) basis and each target qubit in |-|-), 
|— ) basis do not required Hadamard gate which makes 
the procedure for construction NTCNOT-gate simpler. 



B. Motivation and advantages 

Multiqubit quantum controlled phase gate shown in 
Fig. [Uplay a key role in the realization of quantum error 
correction [l^ and implementation of Grover's algorithm 



FIG. 2: (a) Schematic circuit of NTCNOT-gate with qubit 
1 simaltaneously controlling n target qubits. The NTCNOT- 
gate is equavalent to n two-qubit CNOT-gates each have a 
share control qubit 1 but different target qubits 2, 3, or n. 
In this case qubit 1 is deflned in |0) , 1 1) basis while each target 
qubits 2, 3, or n are deflne in |+), | — ) basis, (b) Equivalent 
circuit of NTCNOT-gate in |+), |-) basis. The symbol Z 
represents phase flip of tt on each target qubit. If control 
qubit 1 is in state |— ) than state |— ) at each Z is phase fliped 
as |— ) — > — | — ) while state j+) remain unchanged. However, 
if control qubit 1 is in state |+) , either state |-|-) or state 
|— ) at each Z ramain unchanged. H represents Hadamard 
gate (c) Equivalent circuit of NTCNOT-gate in |0), |1) basis. 
The symbol Z represents phase flip of tt on each target qubit. 
If control qubit 1 is in state |1) than state |1) at each Z is 
phase fliped as |1) — )■ — |1) while state |0) remain unchanged. 
However, if control qubit 1 is in state |0) , either state |0) or 
state |1) at each Z ramain unchanged. It may be noted that 
2(n — 1) Hadamard gates are required in this case. 

for eight objects [l^. Quantum gate with multiple tar- 
get qubits shown in Fig. [5] are of great im por tance in the 
realization of entanglement preparation |14l |. error cor- 
rection ^1, discrete cosine transform jl^, and quantum 
cloning 17]. 

Some interesting scheme for the realization of multi- 
qubit quantum gates have been proposed. In an ear- 
lier studies, J. T Chang et al. [l^ presented three- 
qubit quantum phase gate with four-level atom in a cas- 
cade configuration initially in its ground state through 
a three- mode optical cavity. C. P. Yang et al. pre- 
sented an n-qubit controlled phase gate with supercon- 
ducting quantum-interference devices (SQUIDs) by cou- 
pling them to a superconducting resonator. Recently, 
some schemes are proposed for the realization of a multi- 
qubit phase gate with a fixed phase-shift of tt on each tar- 
get qubit and multiqubit phase gate with random phase- 
shift on each target qubit [2Q-223 . 

Our goal of this work is to realize three qubit controlled 
phase gate shown in Fig. [Ua) and NTCNOT-gate shown 
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in Fig. [2l^a) with four level quantum system in cavity 
or coupled to superconducting resonator. Our proposal 
have several advantages: (i) Decoherence due to sponta- 
neous decay of level |3) is suppressed because the excited 
level 1 3) is unpopulated during the gates operation, (ii) 
The level spacing of the qubit system during the gates 
operations which may cause decoherence is not needed, 
(iii) Operation time for the realization of NTCNOT-gate 
is independent of the number of qubits. (iv) In case of 
flux (SQUID) qubit system each qubit can have much 
longer storage time as discussed in Sec. IV. (v) We do not 
require identical coupling constants of each qubit system 
with cavity mode. Similarly, the detuning of the cavity 
mode with the transition of the relevant levels in every 
target qubit system is not identical, therefore our scheme 
is tolerable to inevitable non-uniformity in device param- 
eters, (vi) Finite second order detuning d — Ac — is 
not required which improves the gate speed by one or- 
der, (vii) Three qubit controlled phase gate shown in 
Fig. [T]is generalized to n-qubit quantum gate with mul- 
tiple control qubits. Intersestly, complexity (number of 
operations) reduces linearly as compare to conventional 
gate decomposition method. In addition, our proposal 
is general and can be applied to various kind four level 
physical system like superconducting devices coupled to 
superconducting resonator and trapped atoms in a cav- 
ity. 



II. SYSTEM DYNAMICS 
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FIG. 3: Desire four level qubit system with four energy level 
jo) , |1) , |2) , and |3). (a) Represents charged qubit system and 
transition freqencies between the levels satisfies the condition 
V2\ > i^io, f32 and ^32 < 1^10. (b) Represents phase qubit sys- 
tem and transition freqencies between the levels satisfies the 
condition satisfies the condition i/io > V21 > ^32- (c) Repre- 
sents flux qubit system and transition freqencies between the 
levels satisfies the condition satisfies the condition P21 > vio, 
1/32 and 1/32 > vio- (d) Represents the SQUIDs qubit system 
and transition freqencies between the levels satisfies the con- 
dition satisfies the condition z/32 < V21 < 1^20 < z^si < ^^30- 
The level |0) and |2) lies in right well of SQUID while level |1) 
lies in left well of SQUID. Their is potential barrier between 
these two wells. 



We consider here four level qubit system which could 
be either natural atoms or artificial atoms as shown in 
Fig. [31 It may be mentioned that Fig. |3] (a) applies 
to superconducting charged qubit [23| . Fig. [31(b) applies 
to phase qubit system [2J, [23 , Fig. [31 (c) applies to flux 
qubit system [23, [2^ and Fig. [31 (d) applies to super- 
conducting quantum interference devices (SQUIDs) [27| . 
Four level structure shown in Fig. [31 (b) could also be 
applied to atoms (2]| . 



A. System-cavity-pulse resonance Raman coupling 

We consider four- level qubit system 1 and 2 coupled to 
single-mode cavity field and driven by classical microwave 
pulse as shown in Fig. [4l (a) and (b). Consider qubit 
system 1 for which cavity mode is coupled to |2)j^ o 
|3)-^ transition but highly decoupled from the transition 
between any other two levels. In addition, microwave 
pulse is applied which is coupled to |l)j^ O transition 
but highly decuple from the transition between any other 
two levels as shown in Fig. 21 (a). Thus, Hamiltonian of 
the system can be written as 
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H = huJca^'a +'^En\n)^{n\ + hgi{a'' \2)^{3\ + H.c.) 

+hn,3{e'^-'\l)^{3\ + H.c), (4) 

where (a) is the photon creation (annihilation) op- 
erator of the cavity mode with frequency ojc, gi is the 
coupling constant of the between the cavity mode and 
\2)^ ■(-> |3)-^ transition of qubit system 1. Rabi frequency 
of pulse is r2i3 and the frequency of the pulse is w^. 

We assume that cavity mode is off-resonance with 
1 2)^ ^ 1 3)^ transition of the qubit system 1 (i.e., Ac = 
LU32 — >> 3i)- Here, Ac is the detuning between 
\2)^ O |3)j^ transition frequency W32 of the qubit sys- 
tem 1 and frequency of cavity field Wc. Microwave pulse 
is off-resonance with ■(-> |3)j^ transition of the qubit 
system 1 (i.e., A^ — W13 — >> rJis). Here, A^ is the 
detuning between o |3)j^ transition frequency wis of 
the qubit system 1 and frequency of pulse oj^. The level 
|3)j^ can be eliminated adiabatically [2g. Thus, for case 
A^ = Ac, effective Hamiltonian in interaction picture 
(assuming h= \) can be written as [131 
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Hr = - 



"13 
Ae 



^«^«|2>i(2| 



(at|2),(l|+i/.c.) 



(5) 



The last two terms describes resonance Raman coupling 
between levels and |2)-^ . For the case flis = gi, initial 
state |2)-^ and |0)^ of the qubit system 1, under 
the Hamiltonian [5] can be written as 

\l)\0)^^e^''[cos{e)\l)\0)^-tstn{d)\2)\l)^], 

|2) ^ e^%cosi9) |2) - zszuiO) |1) |0) J. (6) 

Here, 6 = gft/Ac and |0)^ (|l)c) the vacuum state 
(single photon state) of the cavity field. The state 
|0).|^ |0)c I'smain unchanged under the Hamiltonian [5] For 
the pulse duration ti — TTAc/{2gf) (i.e., 9 = ^) we 
obtained the transformation |0)^ — >■ |2)^ and 
|l)c ~^ |0)c for qubit system 1 and cavity field. 
We will denote this transformation as Gi. 

In the case of qubit system 2, for notation convenience 
we denote ground state (first excited state) as level |1)2 
(10)2) as shown in Fig. |4] (b). Cavity mode is cou- 
pled to 12)2 O 13)2 transition and microwave pulse is 
coupled to 10)2 ^ 13)2 transition of qubit system 2 as 
shown in Fig. |4] (b). In similar fashion, for pulse du- 
ration ^2 = 7''Ac/(2(7|) we obtained the transformation 
10)2 |0), ^ 12)2 |1), and 12)2 ^ 10)2 |0), for qubit sys- 
tem 2 and cavity field. We will denote this transforma- 
tion as G2. However, state |1)2 |0)^ and |1)2 of qubit 
system remain unchanged under transformation G2. 



B. System-cavity off- resonant interaction 
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FIG. 4: (a) System-cavity-pulse resonance Raman coupling 
for qubit system 1. Here, Ac = LJ32 — is the detuning be- 
tween 1 2) J j3)-^ transition frequency u}z2 of the qubit system 
1 and frequency of cavity field uJc, while = ui\s —^^ii is the 
detuning between -H- |3)-^ transition frequency 0^13 of the 
qubit system 1 and frequency of pulse o;^. Both detunings 
are set to be equal (i.e, A^ = Ac) to establish resonance Ra- 
man coupling between level and |2)-^ . Rabi frequency of 
pulse applied is fiia and gi is the coupling constant between 
the cavity mode and |2)j -h- |3)-^ transition of qubit system 1. 
(b) System-cavity-pulse resonance Raman coupling between 
level 10)2 and \2)^ for qubit system 2. Rabi frequency of pulse 
applied is f2o3 and 32 is the coupling constant of the between 
the cavity mode and |2)2 |3)2 transition of qubit system 
2. (c) System cavity off-resonant interaction for qubit system 
k = 2,3, ...,n. Cavity mode is off-resonant with |2)j. -H- |3)^ 
transition of qubit system k with detuning Ac,fe and coupling 
constant gk- 



Consider qubit system fc, for which cavity field in- 
teracts off-resonantly with |2)j, o |3)j, transition (i.e., 
Ac,/c = ujc — W32 >> gk) while decoupled from any tran- 
sition between the other levels as shown in Fig. |4] (c). 
Here, Ac,fc is the detuning between |2)j, o |3)j, transi- 
tion frequency 1032 of qubit system k and frequency of 
cavity field uJc while gk is the coupling constant between 
resonator mode and |2)^, o |3)j. transition. The effec- 
tive Hamiltonian of system in interaction picture can be 
written as [2^ 

^i = ^(|3).(3|-|2),(2|)ata. (7) 
Ac,fe 

In the presence of single photon in cavity, the evolution 
of initial state |2) |I)^ and |3) is given by 

|2)Jl).^e^^^/^-|2),|l)^, 

|3)Jl),^e-^^*/'^-M3)Jl),. (8) 

It is clear that phase shift of e*^^*/'^^-'' (e^^s'*/'^-.'') is 
induced to the state |2)^ |1)^ (|3) 1^1)^) of qubit system k. 
However, states \2)f, |0)^ and |3)j^, |0)^ remain unchanged. 



C. System-pulse resonant interaction 

Assuming that microwave pulse resonant to |j) — > |2) 
transition of each qubit system is applied. Here, j = 1 
for qubit system 1 and fc, while j = for qubit system 2. 
Then, the evaluation of state is given by [19| 

\j) ^ co.s{n,2T) \j) - te-'^sin{n,2T) |2) , 

|2) ^ cos{nj2T) |2) - ie^'^sm(%2r) |j) , (9) 

where Qj2 is the Rabi frequency between two levels 
\j) and |2). Where as r represents interaction time 
of qubit system with microwave pulse and (p is the 
phase associated with microwave pulse. For pulse du- 
ration T = 7r/(2f2j2) and phase f — 7r/2, transformation 
|2) (b)) — > |j) (— |2)) is obtained which is denoted by R. 
For the case of phase ip = —tt/2, we obtained the trans- 
formation |2) (|j)) -)> - |j) (|2)) denoted by R'' . It may be 
mentioned that resonant interaction of microwave pulse 
with qubit system can be carried out in a very short time 
by increasing the Rabi frequency of pulse. 
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III. IMPLEMENTATION OF MULTIQUBIT 
CONTROLLED GATE 



The goal of this section to demonstrate how three- 
qubit quantum phase gate and NTCNOT-gate can be 
reahzed based on system dynamics described in Sec. II. 



A. Three-qubit controlled phase gate 

We consider qubit system 1, 2 and k (with A: = 3) for 
the implementation of three-qubit controlled phase gate 
shown in Fig. 3]. For each qubit system two lowest energy 
levels |0) and |1) represents logical state of each qubit 
while other higher energy levels |2) and |3) are utilized 
for gate realization. We assume that cavity is initially in 
vacuum state |0)^. Three-qubit controlled phase gate is 
realized in the following steps 

Step (i): Apply the transformation Gi to qubit system 

1 for time ti. Namely, when qubit 1 is initially in state 
|l)j^, a photon is emitted in cavity. However, the state 
|0)-^ |0)^ remain unchanged under transformation Gi. 

Step (a): Apply transformation R to qubit system 1 
and i?^ to qubit system 2, simultaneously. In this step we 
set T = tt/ (2rio2) = 7r/(2r2i2) by adjusting the intensities 
of two pulses. 

Step (iii): After the above operations level \2)^ of qubit 
system 1 is unpopulated. The level |0)2 of qubit system 

2 transform to level 12)^. Apply transformation G2 (for 
time duration ^2) to qubit system 2 which will absorb the 
single photon from the cavity. However, if qubit system 
2 is in state |1) single photon from the cavity will not be 
absorbed. 

Step (iv): Apply transformation i?^ (for time dura- 
tion r) to qubit system fc = 3. After this operation, 
when cavity is in single photon state level |2) of both 
qubit system 1 and 2 are unpopulated. Under this con- 
dition cavity field interacts off- resonantly to 12)3 — > 13)3 
transition of qubit system 3. It is clear from Eq. ^ 
that for tk — {nAc,3)/gl the state [2)3 |1)^ of qubit 
system 3 changes to — [2)3 However, the states 

10)3 |0)^, 10)3 and 12)3 |0)^ of qubit system 3 remain 
unchanged. Finally, apply transformation R (for time 
duration r) to qubit system 3. 

Step (v): Apply transformation G2 (for time duration 
^2) to qubit system 2. 

Step (vi): Apply transformation R^ to qubit system 1 
and R to qubit system 2, simultaneously for time dura- 
tion T. 

Step (vii): Apply the transformation Gi to qubit sys- 
tem 1 for time ti. Asa result, qubit 1 is transformed 
back to state and cavity field returns to its original 
vacuum state. 

Above all operations are schematically presented in 
Fig. [S] The states of the whole system after the above 
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FIG. 5: Scematic diagram for the implementation of three- 
qubit controlled phase gate. Gi represents the system-cavity- 
pulse resonance Raman coupling between level and |2)j 
for qubit system 1 while G2 represents system-cavity-pulse 
resonance Raman coupling between level |0)2 and \2)^ for 
qubit system 2. Gtt is system cavity ofT-resonant interaction 
where as R and represents system pulse resonant interac- 
tion. 



operations are summarized below 



1 100) 
jlOl) 
jllO) 
|111) 
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1111) 
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1111) 11) 
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1 200) 
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1210) |i: 

-1211)11 



1100) |1), 

1101) |1), 
1110) |1), 

■|111)|1). 



(10) 



Here, state \abc) is the abbreviation for the states \a)-^, 
\b)^ and \c)^ of qubit (1, 2, and 3) with a, 6, c G [0,1,2]. 

On other hand, states |000) |0)^, |001) |0)^, |010) |0)^, 
and |011) |0)^ remain unchanged. This is because state 
|0)j of the qubit system 1 is not effected by the applica- 
tion of transformation Gi that is no photon in cavity is 
emitted when qubit 1 is in state |0)j^ . Hence, it is clear 
from Eq. pH)) that three-qubit controlled phase gate was 
achieved with three qubits (i.e., control qubit 1 and 2, as 
well as target qubit 3) after the above operations. 

Present proposal provide simple way to realize the Tof- 
foli gate shown in Fig. [Ijb) . It is well known that at least 
six two qubit controUed-NOT-gates and ten single-qubit 
gates (i.e., two Hadamard, one phase , and seven 7r/8 
gates) are required to construct a Toffoli gate by conven- 
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tional gates decomposition methods [3C|. The two qubit 
CNOT-gate is equivalent to two Hadamard gate and sin- 
gle two-qubit phase gate. If we assume that realization 
of single-qubit gate and two-qubit phase gate requires 
only one step operation, than, using conventional gate 
decomposition method, at least 28 steps will be required 
to realize Toffoli gate. However, present proposal requires 
only 9 steps i.e., 7 steps for three-qubit phase gate plus 
two steps operations for two Hadamard gate. 

Our scheme can easily be generalized to n-qubit con- 
trolled phase gate with multiple control qubits. For this 
purpose, (i) Apply transformation Gi and R to qubit 1. 
(ii) Apply transformation i?^ and G2 to qubit system 
2,3,...,n — 1. (iii) Apply transformation i?^, G^, and R 
to last qubit. (iv) Apply transformation G2 and R to 
qubit system 2,3, ...,n— 1. (v) Apply transformation R^ 
and Gi to qubit 1. Hence, n-qubit controlled phase gate 
can be achieved by a sequence of operations which are 
summarized as 



i=2 



C/^ = Gi (g) i?'!' (g) J]^ [i? (K) 62] (8)(i?(8)G^ 

Tl-l 

n— 1 . j 

(gi?^)® ]^[G2«)i?^] ®i?«)Gi, (11) 



where Hn-^i G ~ G 



i=2 
3 

5 G ( 



n—1 - 2 

G while nr=2 G = 



G (g) ... g) G (g) G. Total number of steps for n-qubit 
quantum gate with multiple control qubits will be 4n — 5. 
Assuming that realization of basic gates requires only one 
steps operation, than, using conventional gate decompo- 
sition method, required number of steps will be 22n — 75 
(30| . Hence, our scheme reduces the number of steps 
(complexity) linearly as for as number n of the qubit is 
concerned, as compare to conventional gate decomposi- 
tion method. 



B. NTCNOT-gate 

In order to implement NTCNOT-gate we consider 
qubit system 1 (as shown in Fig. 2] (a)) initially pre- 
pared in state (|0);^ -I- |l)j) /^/2. In this case we consider 
n—1 qubit system of type k shown in Fig. |4] (c) with 
k = 2, 3, ...n . Each qubit system k is initially prepared 
in state |0)j.. In the new rotated basis for qubit system 
k, the state of the whole system can be written as 



^(|0)i + |l)i)' 



11(1+)^ + 1-)^^)' 



|1)-^, a photon is emitted in cavity. However, the state 
|0)-^ \0)^ remain unchanged under transformation Gi. 

Step (ii): Apply transformation R to qubit system 1 
and i?^ to each qubit system k for time duration t, si- 
multaneously. As a result transformation 
|a)j, (|fe)j.) is obtained for each qubit system k. Here, 
|a), = l/y2(|0), + |2),) and |6), = l/\/2(|0), ~ |2),).^ 

Step (iii): After above operations, when cavity is in 
single photon state level \2)^ and level |3)-^ of qubit 
system 1 is unpopulated. Under this condition cavity 
field interacts off- resonantly to |2)j, — >■ |3)j. transition of 
each qubit system k. It is clear from Eq. (|5]) that for 
tk — (7rAc,fc)/f/^ the state |2)^ |1)^ of each qubit system 
k changes to — |2) |1)^. That is in the presence of single 
photon in cavity the state |a)j, |1)^ of each qubit system 
k changes to 1^)^ |l)j, while |6)j, |1)^ of each qubit system 
k changes to |a)j. |1)^. However, the states \a)^ |0)^ and 
1^)2 \^)c remained unchanged. 

Step (iv): Apply transformation R^ to qubit system 
1 and R to each qubit system k for time duration r, 
simultaneously. 

Step (v): Apply the transformation Gi to qubit system 
1 for time ti.. As a result, qubit 1 is transformed back to 
state |l)j^ and cavity field returns to its original vacuum 
state. 

After the above operations, one can easily see that 
controUed-NOT gate of one qubit simultaneously con- 
trolling n qubits described by Eq. ^ and Eq. ^ was 
achieved with n qubit system (i.e., control qubit 1 and 
target qubit systems k = 2,3, ...n ). 

In order to understand better realization we consider 
example of three-qubit case. It can be checked that states 
of the whole system after the above operations are sum- 
marized below 



|1 + +>|0), 
|1 + ->|0), 

|i-+>|0), 
|i-->|0), 

|l^'^'>|l)c 

|la&)|l). ^ 
\laa) |1), 



|2 

|2-- 
|2- 
|2--) 
|2-+) 
I2 + -) 
12 -I- -f) 



+)|1> 
->|1> 

->|1) 

|l>c 
|l>c 

Il)„ 



2 



|laa>|l>c 

|la^>|l)c 3 
\Ua) |1), ^ 

|i-->|0), 
|i-+)|o), 

|1 + ->|0), 
|1 + +>|0),. 



(13) 



(12) 



Hence, it can be concluded from Eq. (|13p that three- 
qubit controUed-NOT gate of one qubit simultaneously 
controlling 2 qubits with fc = 2, 3 was achieved with 3 
qubit system (i.e., control qubit 1 and two target qubit 
systems 2 and 3). For fc = 2 in Eq. ^ and Eq. ^ our 
scheme reduces to two-qubit controUed-NOT gate which 
can be used to implement two qubit Deutsch-Jozsa algo- 
rithm as described below 



fe=2 



where, |±)j, = l/'\/2(|0)^± |1)^). Operations for realizing 
NTCNOT-gate are described as fallow. 

Step (i): Apply the transformation Gi to qubit system 
1 for time ti. Namely, when qubit 1 is initially in state 



1. DEUTSCH-JOZSA ALGORITHM 

Deutsch-Jozsa algorithm is designed to distinguished 
between the constant and balanced functions on 2" in- 
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puts [3l|. For constant function, the function f{x) ~ 
constant for all 2" inputs. For the balanced function, 
the function f{x) — for half of all possible inputs, and 
f{x) = 1 for other half. A classical algorithm will need 
2"'/2-|-l queries to determine whether function is constant 
or balanced since there may be 2"/2 O's before finally a 1 
appears, showing that function is balanced. In contrast, 
the Deutsch-Jozsa algorithm requires only one query. 

Next, we discuss the scheme to implement two-qubit 
Deutsch-Jozsa algorithm using four-level qubit system 
shown in Fig. [3] coupled to cavity or resonator. The 
qubit system 1 shown in Fig. |4l[a) represents query qubit 
while qubit system k = 2 shown in Fig. |31Jc) represents 
auxiliary qubit. We prepare the two qubit system in the 
state IV-) = 1/V2{\0)^ + |1) J |1)2 which can be written 
in rotating basis for qubit system fc = 2 as 



->2), 



(14) 



The function f{x) is characterized by the unitary map- 
ping transformation Uf, and \x,y) — >■ \x,y (B f(x)), where 
® is a symbol of addition modulo 2. After unitary trans- 
formation [//, initial state of the system changes to 

|0), + (-1)^(1) |1) J ® (1+)^ - l-)^). (15) 

There are four possible transformations: (i) Uf^i corre- 
sponding to /(O) = /(I) = 0; (ii) Uf^2 corresponding to 
/(O) = /(I) = 1; (iii) [//,3 corresponding to /(O) = 
and /(I) = 1; (iv) [//_4 corresponding to /(O) = 1 
and /(I) — 0. Then Hadamard gate applied on query 
qubit. As a result, the state of query qubit becomes 
1/(0) ® /(I)). If f{x) is constant, the state of query qubit 
becomes |0)-^. On other hand, if f{x) is balanced, the 
state of the query qubit becomes Therefore, a mea- 
surement on query qubit provides the desire information 
whether the function f{x) is constant or balanced. The 
Uf^n operations are applied to the state jV') as fallow 

Uf^i operation. This is an identity operation. Both 
qubit system are kept far off with the cavity field and 
microwave puis. As a result system remain in the state 

IV'). 

U f^2 operation. We first apply two-qubit controlled 
NOT-gate as described above. Next, we apply single- 
qubit rotation |0) |1) and |1) — — |0) on qubit system 
1. Then we repeat two-qubit controlled NOT operation 
and perform the single-qubit rotation |0) — > — |1) and 
|1) — >■ |0) on qubit system 1. Finally, we obtained 



IV')2=o(-|0)i-|l)i)®(|+)2-|-)2)- 



(16) 



U operation. We apply two-qubit controlled-NOT 
operation. Then state of system evolved to 

IV')3 = ^(|0)i-|l)i)®(|+)2-|->2)- (17) 



f//,4 operation. We apply single-qubit rotation |0) — 
|1) and |1) — )■ — |0) on qubit system 1. Then we perform 
controlled-NOT operation. Finally, we apply single-qubit 
rotation |0) — >■ — |1) and |1) — > |0) on qubit system 1. 
The resultant state becomes 



1 



(-|0)i + |l)i)®(|+)2-|->2)- 



(18) 



In this way we obtained the unitary mapping transfor- 
mation Uf of Eq. ()15|) . After Hadamard transformation 
on qubit system 1, if state of qubit system 1 becomes 
|0)j^, the function f{x) is constant. On other hand, if 
state of qubit system 1 becomes the function f{x) 
is balanced. 



IV. POSSIBLE EXPERIMENTAL 
IMPLEMENTATIONS 

In this section, we give a discussion on experimen- 
tal possibilities of three-qubit controlled phase gate and 
NTCNOT-gate described in Sec. III. Total operation 
time of three-qubit controlled phase gate is given by 



rgcp = 2ii + 2*2 + 4- + 4r 



,7rA 



Secondly, total operation time for NTCNOT-gate is given 
by 



= 2ti + 2t + tk = 2(^) + 2( 



2gi' ' -^2f^i2^^^ gl ^ 
(20) 

The operation time T^p and Tntcnot should be shorter 
than (i) energy relaxation time of level |2) (it may 
be mentioned that level |3) is unpopulated during the en- 
tire operations) , and (ii) the life time of the cavity mode 
— Q/2'!Th'c, where, Q is loaded quality factor of the 
cavity and i>c is the resonator frequency. In principle, 
these requirements can be achieved through (i) reduc- 
ing operation time by increasing coupling constants and 
Rabi frequencies, (ii) increasing by employing high 
Q cavity or resonator, and (iii) choosing qubit system 
(e.g., atoms) or designing qubits (e.g., superconducting 
devices) such that the energy relaxation time of level 
1 2) is sufficiently long. 

Here, we consider without loss of generality (?i ~ 52 ~ 
gk ~ g. On choosing Ac ~ Ac, 3 ~ Ac,fc ^ lOg, and 
fli2 ~ 10(7, total operation time required for the gates 
implementation would be r3cp ~ 30.27r/(7 and Tntcnot ~ 
20TT/g. For rough estimate, we assume g/ir ~ AAOMHz, 
which could be achieved for superconducting qubits cou- 
pled to a one-dimensional standing-wave coplanar wave 
guide (CPW) transmission resonator [s^. For g cho- 
sen here, we have r3cp ~ 0.068^s and Tntmot ~ 0.045/is, 
much shorter than ~ l^s (33| . In addition, consider 
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FIG. 6: An rf SQUID with first four energy levels. Magnetic 
dipole coupling between two ground levels |0) and |1) is much 
smaller than that between any other levies due to potential 
barrier between two wells. Transition frequencies between the 
excited levels and ground levels are 1/30 ~ 24:.4GHz, 1/31 ~ 
21AGHz, U12 ~ 16.5GHz, U20 ~ IQ.^GHz, which are much 
larger than transition frequency 1^32 ~ A.9GHz. 

resonator with frequency Vc ^ 3GHz and Q ~ 10^ 8], 
we have ~ 5.3/is which is much longer than oper- 
ation times Tcp and Tntcnot- It may be mentioned that 
superconducting coplanar wave guide resonator with a 
(loaded) quality factor Q ~ 10^ have been experimen- 
tally demonstrated j34i] . 

As a quantitative example of our technique, consider a 
superconducting quantum interference devices (SQUIDs) 
as qubit system. For SQUID, level structure is straight 
forward to implement by changing external control pa- 
rameters e.g., magnetic flux (f>x '27i|. Fore example, con- 
sider rf SQUID shown in Fig. [5] with junction capac- 
itance C = 90/F, loop inductance L = lOOpH, junc- 
tion's damping resistance R ^ IGO, potential shape pa- 
rameter = 1.12, and external flux (j)^ = 0.4995(/)o. 
Here, (f>Q = h/2e is flux quantum. It may be men- 
tioned that SQUIDs with these parameters are avail- 
able currently [23|. With these choices, decay time of 
level |2) would be ~ lOO/iS, the |2) -J> |3) cou- 
pling matrix element is (f>32 ~ 7.8 x 10"^, and |2) — > 
|3) transition frequency is 1^32 ~ A.QGHz. We choose 
cavity mode frequency Vc = Wc/(27r) — S.GGHz and 
Q ~ 10^ we have k~i ~ 4.42^s. The SQUID-cavity 
coupling constant for |2) — > |3) transition is given by 
g = {l/L)y/ujc/'2fiQh(f)32(f)a JgBc{r).dS. Here, S is the 
surface bounded by the SQUID ring and Bc(r) is the 
magnetic component of cavity mode in the SQUID loop. 
For standing wave cavity, Bc{z) — [iQ\plfV cos kz, where 
fc, V , and z are wave number, cavity volume, and the cav- 
ity axis, respectively. For 10 x 1 x Xmrr? standing-wave 
cavity and a SQUID with 40 x AOfirri^ loop located at 



antinode of the cavity mode. Simple calculation shows 
that g ~ 4.3 x lO^s"^ for which time required for (i) 
three-qubit phase gate would be r^cp ~ 0.219/is and (ii) 
for NTCNOT-gate would be T,uc,iot ~ 0.146^s. These 
implementation times are much shorter than ^-iid 
K^^. Moreover, we have an additional advantage in case 
of flux-qubit system that is tunneling between the lev- 
els |1) and |0) is not needed during the gates operation. 
Therefore, potential barrier between levels |1) and |0) can 
be adjusted priory such that decay from level |1) becomes 
negligibly small [23j[lE|- As a result, each qubit can have 
much longer storage time. 

During the operation of step (ii) or (iv) or (vi) 
for three-qubit phase gate and of step (ii) or (iv) for 
NTCNOT-gate, a single photon is populated in the cavity 
mode and state |2) of each qubit system is occupied. The 
unwanted system-cavity-pulse resonance Raman coupling 
and system-cavity off-resonant interaction between res- 
onator mode and |2) — > |3) transition of qubit system in- 
duces as accumulated phase shift to state |2) of each qubit 
system, which will effect the desire gate performance. 
However, when r << ti,t2,tk this unwanted phase shift 
is sufficiently small and can be neglected. Note that for 
ili2 = r2o2, we have t = 7r/(217i2), ti = 7rAc/(2gJ), 
t2 = 7rAc/(2(7|), and tk — TrAck/gl- Thus condi- 
tion turns into fti2 >> 2g1/ Ac,2g2/ g^/ ^c.k which 
can be met by increasing the Rabi frequency of pulse 
(i.e, by increasing the intensity of resonant pulse). For 
Ac,fe = lOgfe, the occupation probability of level |3) of 
target qubits is approximately 0.04 which reduces the 
gate error [27|. 



V. CONCLUSION 

Before conclusion, we should mention some interesting 
previous proposals on multiqubit gates [T8l - [2l| . However, 
the present alternative approach is different from these 
previous proposals. The proposal [l^ employs three- 
mode cavity. The present scheme requires only single 
mode cavity or resonator which makes the technical re- 
quirements comparatively easy. As compare to the pro- 
posal [1^, decoherence due to adjustment of level spac- 
ing (which is also undesirable in experiment) is avoided. 
In our previous proposal [13] level ]3) is populated dur- 
ing the gate operations, the spontaneous decay from this 
level may cause decoherence. In the present scheme, level 
]3) is unpopulated during the entire operations, hence 
suppressed the spontaneous decay of this level. The pre- 
vious proposal t21||] requires 8 steps operations to im- 
plement NTCNOT-gate, assuming that Hadamard gate 
can be implemented in single steps and also required fi- 
nite second order detuning S — Ac — A^. However, us- 
ing present proposal NTCNOT-gate can be implemented 
with 5 steps without Hadamard gate and finite second 
order detuning is not required which improves the gate 
speed by one order. 

We have proposed a scheme for realizing three-qubit 
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controlled phase gate and NTCNOT-gate as shown in 
Fig. dD) and Fig. ^ with three types of interac- 
tions. These interactions are system-cavity-pulse reso- 
nance Raman coupling, system-cavity off-resonant inter- 
action, and system-pulse resonant interaction. Our pro- 
posal can be applied to various kind of physical system 
with four-level configuration. For different systems, fre- 
quency regimes of cavity mode could be different, e.g., 
optical cavities in case of atoms and microwave cavities 
in case of superconducting qubits. However working of 
these gate in light of experimental errors requires a rather 
lengthy and complex analysis which should be further in- 



vestigated. 

We have shown that our proposal have several advan- 
tages: (i) Decoherence due to spontaneous decay of level 
1 3) is suppressed because the excited level |3) is unpopu- 
lated during the gates operation, (ii) The level spacing of 
the qubit system during the gates operations which may 
cause decoherence, is not needed, (iii) For the quantum 
gate with multiple control qubit, the number of steps 
(complexity) reduces linearly for number n of the qubit, 
as compare to conventional gate decomposition method, 
(iv) Operation time for the realization of NTCNOT-gate 
is independent of the number of qubits. 
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